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ABSTRACT 


The  determination  of  all  the  extreme  points  of  a  given  convex 
polyhedron  P  —  R  n  generally  requires  a  substantial  amount  of 
computations;  this  note  presents  a  conceptually  simple  algorithm 
for  this  purpose.  Unlike  other  methods,  the  procedure  generates 
only  those  basic  solutions  which  are  extreme  points  (i.e.,  only 
feasible  basic  solutions). 

More  generally,  this  approach  is  able  to  generate  all  the 
faces  of  any  dimension  k  (0  <  k  <  n)  ,  that  is  all  those 
k-dimensional  subpolyhedra  which  lie  on  the  boundary  of  the  given 


polyhedron  P  . 


Introduction 


Polyhedral  sets  are  the  most  widely  used  constrained  sets  in  mathe¬ 
matical  programming  and  they  are  usually  defined  by  a  system  of  linear 
(in)equalities.  In  fact,  other  characterizations  of  these  sets  are  often 

impractical  in  the  convex  programming  context. 

However,  as  mathematical  programming  penetrates  into  the  darker  and 

less  structured  areas  of  non-convex  programs  (concave  programming  or  dis¬ 
crete  (in  particular,  integer  or  zero-one)  programming,  for  instance),  more 
information  on  the  polyhedral  set  P  of  feasible  solutions  is  required. 

A  classical  question  is  that  of  finding  the  vertices  of  P  ,  and  seems 
very  difficult  to  answer  practically  in  large  dimensional  vector  spaces. 

One  may  also  be  interested  in  the  complete  face  structure  of  P  ,  that  is, 
in  finding  a  characterization  of  each  k-dimensional  face  F  of  P  ,  for 
0  <  k  <  n  .  (Clearly  the  quest  for  vertices  is  a  special  case  of  the  latter, 
since  they  are  O-dimensional  faces  of  P.) 

This  paper  presents  a  simple  approach  to  the  determination  of  the  face 
structure  of  a  polyhedron.  An  algorithm  is  presented,  which  generates  all 
the  desired  information  concerning  the  complete  lattice  of  faces  of  a  poly¬ 
hedron  P  ,  in  the  form  of  a  non-redundant  facial  arborescence. 

Somj  applications  are  mentioned  and  described  in  general  terms  in  the 
last  section.  One  particular  case  is  that  of  general  quadratic  programming, 
which  is  the  object  of  the  follow-up  paper  [4], 


Section  1.  Minimal  Sets. 

1.1  Consider  the  polyhedral  set  P  ,  defined  by  the  following  system  of 
inequalities 


x.  =  a.  + 
1  10 


E  a.  . 
jeN  1J 


x.  >  0 

J  ~ 


i  e  M 


(1) 


with  non-basic  index  set  N  —  {l ,2 , . . .  ,n}  C  M  =  (l ,2 , . . . ,n ,n+l , . . . ,n+m} 


-  The  matrix  A  contains  an  n-by-n  identity  submatrix,  which 
correspond Ss  to  the  constraints  x^>0,  jeN-~M 

-  For  simplicity,  we  only  consider  here  the  case  where  P 
has  full  dimension  n  and  is  bounded. 

Def inition:  A  subset  1 C  M  is  called  minimal  if 

ixj  x  ^  >  0 ,  Vi  e  l}C^{x|x^>0,  Vi  e  m}  =  P 

and  for  every  iQ  e  I  ,  there  exists  a  point  x  such  that 

x.  <0  and  x.  >  0  ,  Vi  e  I  -  f i  ]  (2) 

1  i  —  k  o 

o 

Property  1:  The  constraints  >  0  ,  Vi  e  (M  -  I)  are  redundant 
and  one  has 


O) 


P-  (x  x.  >0  ,  Vi  e  l] 

L  — 


2 


Property 

Proofs : 

to  1 

to  2 


For  every  i  e  I  ,  the  hyperplane 
(n-1) -dimensional  facet  of  P  . 


x.  =0  contains  a 

l 

o 


By  definition  one  has 

V  q  e  (M-l  ):  >  0  ,  V  x  e  (x  j  x.  >  0  ,  V  i  6  I  }  q.e.d. 


The  set 

D  = 

lx  1  x.>0,VieI 

> 

-  {iQ}  J  contains  a 

point 

x  with 

x . 
l 

o 

<  0  ,  because  i  el 
o 

.  Furthermore  there 

exists 

a  point 

X 

in  the  interior  of  D 

with  x.  =  0  ;  such  a 
xo 

point 

may  be  constructed  by  choosing  an  arbitrary  interior  point 

x  e  Int(P)  and  intersecting  the  line  x  -  x  with  the  plane 

x.  =  0  ,  yielding  x  .  Since  P  c  D  one  has  x  e  Int(D) 

1o 

due  to  the  convexity  of  D  ,  i.e.  x  =  (1-p)  x  +  p  x  ,  with 
0  <  p  <  1  .  Moreover,  since  PC  DC  K n  ,  the  polyhedral  set 
D  has  the  same  dimension  n  as  the  set  P  . 

Take  now  a  (small  enough)  open  n-dimensional  ball  B(x)  ~  Int(D)  , 

containing  x  ,  and  consider  the  intersection 

B  =  fx  e  B(x)  |  x.  =  0  }  :  by  construction  b  is  a 

o 

(n-1) -dimensional  open  ball  C  D  ,  and  it  lies  in  the  inter¬ 
section  of  P  with  the  hyperplane  x^  =  0  ,  i.e.  in  a  facet  of 

o 

P  ,  which  has  dimension  (n-1)  .  q.e.d. 
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1.2  The  concept  of  minimal  set  of  inequalities  I  provides  for  the 
basic  tool  of  an  algorithm  for  the  face  decomposition  of  P  ;  I 
can  be  obtained  by  the  following: 

PROCEDURE  MIN  SET 

Give  a  system  of  inequalities  (1),  i.e.,  a  matrix  A  with  (rrt-m)  rows 
and  (n+1)  columns,  containing  an  identity  submatrix;  the  procedure  MIN SET 
determines  the  minimal  set  I M.  A  primal  feasible  linear  programming 
tableau  for  the  system  (1)  is  required  to  start  MIN  SET.  During  the  exe¬ 
cution  of  MIN SET,  the  elements  i  e  M  are  selected  ore  after  the  other, 

o 

and  the  corresponding  row 


x .  =  a .  +  E  a .  .x . 

i  1  0  i  j  l 

o  o  jeN  o 


is  (momentarily)  chosen  as  objective  function.  Optimization  of  the  following  L.P. 

Minimize  x.  ,  subject  to  x .  >  0  ,  Vi  e  R  -  f  i  furnishes  a  minimal 
l  i  —  k  o 

o 

value  x.  ;  the  index  set  R  is  determined  by  the  procedure  and  satisfies 


I  l  RCM.  If  x.  <0  then  one  has  l  s  I  ,  by  definition  of  the  minimal 
1  o 

o 

set  I  ;  if  x .  >0,  then  the  constraint  x  >  0  is  considered  redundant 

i  l  ~ 

o  o 

and  is  disposed  of;  the  cast  x.  =  0  is  treated  separately. 

o 


MINSET :  1  Set  R  =  M  and  I  =  I  =0 

o 

2  Choose  a  basic  index  r  c  ^R  -  (R  H  N )  j  ; 

-  if  R  =  0  then  STOP! 

-  if  R  ^  0  but  R  c  N  ,  then  change  the  basis  (and 
the  non-basic  set  N  )  by  choosing  a  positive  pivot 
which  preserves  primal  feasibility;  if  there  exists 
no  such  element  in  the  current  tableau  then  replace  1 
by  I  (J  R  and  S10P! 

Set  R  =  R  -  {r}  and  consider  the  r1"*1 


x 

r 


a 

ro 


+ 


! 

jeN 


a 

r  ) 


3 


row 


4 


4  Solve  by  L.P.  optimization  the  problem 

Minimize  x^  ,  s.t.  >  0  ,  Vie(RL)I) 

5  If  the  minimal  value  of  x  is  >  0  then  go  to  2. 

r 

If  the  minimal  value  of  x  =0,  then  replace  I  by  I  U  [  i  } 

r  o  o  o 

and  go  to  2. 

6  Replace  the  I  by  I  U  {r}  and  go  to  2. 


Because  the  procedure  MINSET  will  be  frequently  called  in  the  algorithm 
it  is  well  worth  noting  the  following  remarks  to  speed  up  its  execution: 

1)  In  order  to  minimize  the  number  of  pivotal  operations,  the  choice 

of  r  in  Step  2  should  correspond  to  a  row  x^  with  the  smallest 

possible  number  of  negative  elements  a  (usually  with  just 

one  a  .  <  0  )  . 
rj 

2)  Suppose  that,  in  the  course  of  the  (primal)  optimization  of 

Step  4,  a  row  x  ,  s  s  R.with  a  >  0,  Vi  is  found  in  the  current 
-  s  sj  ~  J 

basic  sec;  then  s  may  be  immediately  discarded  from  the  set  R. 

Proof:  In  the  current  basis,  the  condition  a  ,  >  0  is  the  optimality 

-  sj  - 

criterion  of  the  L.P. 

Minimize  x  ,  subject  to  x.  >  0,  Vi  e  ((R  U  1)  ■  t  s } )  . 

3  1 

Thus,  one  may  bypass  the  minimization  of  x^  and  go  to  the  steps 

5  and  6,  with  x  =  a 

s  so 

3)  Similarly  for  the  column  of  a  non-basic  variable  with  j  e  R, 
one  may  eliminate  the  element  j  from  the  set  R.  Define 


KG 

min  >,  -  i  v.  -  basic  index  (also  k  =  r) 

k/a.  .>  0  akj 
kj 


Then  j  may  be  eliminated  when  A  >  0. 


(A  particular  such  case  is  wnen  <  0,  Vk) . 

Proof :  The  preceeding  condition  is  such  that  the  solution  remains  feasible 


when  Xj  assumes  a  negative  value  0  >  x^  >  -A. 


Let  us  now  show  that  the  procedure  MIN  SET  does  indeed  determine  a 

minimal  set  X  .  : 

min 

(i)  By  construction,  every  point  x  c  P1  *  {x  (  >  0,  i  e  1^  } 

satisfies  x.  >0,  Vi  e  (M-I  ); 

i  —  min 

hence,  P'c  P  =■  £x  |  x^  >  0,  Vi  e  M}  . 

(ii)  For  every  i  e  I  .  ,  the  procedure  MINSET  constructs  a 

o  rain 

point  x  with  x.^  <  0  and  >  0  ,  Vi  s  (R  u  I)  where  R  and  I 
o 

are  the  current  sets  of  seep  4;  since  (RUE)  3  by 
construction,  one  has  x .  >  0,  Vi  e  I  .  -  f i  ] . 

(iii)  The  set  Iq  (4  0)  indicates  that  the  given  system  (1)  is 

degenerate,  i.e.,  that  some  k-diraensional  faces  (0  <  k  <  n) 

of  the  polyhedron  P  are  contained  in  more  than  (n-k) 

hvperplanes  x^  =  0,  i  e  M;  this  situation  does  not  affect 

the  minimal  property  of  *mtn»  but  the  identification  of  the 

elements  of  I  is  important  in  order  to  eliminate  redun- 
o 

dancy  of  the  facial  decomposition.  (See  section  3.5). 

(iv)  Property  3:  In  the  non-degenerate  case  (Iq  *  0)?  the  mini¬ 
mal  set  is  unique. 


Proof :  For  every  i  e  (M-I)  one  has  x  >  c  >  0,  Vx  c  P  since  Iq  *  0 
by  hypothesis;  this  is  a  well  defined  criterion  which  divides 
M  uniquely  in  two  disjoint  subsets  (M-l)  and  I. 
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Section  2:  The  Faces  of  P  . 

At  the  beginning,  one  applies  MINSET  to  the  system  (1),  that  is,  to 
the  given  polyhedron  P  in  order  to  determine  the  set  I  (degeneracy  will 
be  considered  separately  in  3.3).  But  from  the  theory  of  polyhedral 
sets,  one  knows  that  every  face  F  of  the  polyhedron  P  is  a  polyhedron. 

Thus,  MINSET  can  be  applied  to  the  faces  F  of  P  as  well.  In  particular, 
to  the  faces  Fd^),  Fd^.i^)  , -  where 

F(i  ')  =  fx  e  P  I  x.  —  0  ,  i  c  I  }  ,  with  a  corresponding  minimal  set  TO.,)  c-  I  , 
l  ■  1  i,  i  o  1 


F(i.,i^>  =  {x  e  P 


il  e  I,  i2  e  1(1,)  c  i  }  "  I(i1,i2) 


I(l1)  , 


etc 


A  sequence  of  faces  F(i^),  F(i^,i2),.. 
with  the  following  properties. 

Property  3 :  For  k  +  s  <  n  ,  one  has 

Id.  ,i- » . . .  ,i  )  I(x. ,x~ , 
1  l  s  i  l 

and  F(i  ,i„,...,i  )—  F(i  ,i  , 

I  C.  3  1  4 


,F(i1,...,i<?)  is  generated 

"’Wl . in-k) 

*  *  * is ,J  s+1 ’ ’ *  * ’ in-k)  * 


Proof: 


by  construction. 


In  conclusion,  one  sees  that  repeated  use  of  the  procedure  MINSET, 

leads  to  the  construction  of  an  aborescence  with  initial  node  ?  (n-dimen- 

sional  face)  itself;  at  the  level  below,  one  finds  all  the  (n-l)-diroen- 

sional  faces  of  P  (one  for  each  i  s  I  );  then,  below,  the  faces  of  these 

faces  (i.e.,  the  (n-2)-dimensional  faces  of  P)  etc.,...  At  every  node 

(i.e.,  face  F(i,,...,i  ))  the  minimal  set  )  determines  the 

1  q  1  q 
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branches  (how  many  and  which)  leading  Co  Che  Level  below.  At  the  lowest 
level,  one  ultimately  finds  the  vertices  of  P. 

The  next  section  presents  an  algorithmic  construction  of  this  arbor - 
escence. 


i 


8. 


Section  3:  An  exhaustive  arborescence  for  the  faces  of  P  . 

3.1  The  following  procedures  TREE,  FACE2D,  S1MFACE  and  BACKTRACK  generate 
one  by  one  a  list  of  n-arrays  called  VERTEX  =  {i  ,i^ , . . .  ,i^}  . 

In  this  version,  the  algorithm  requires  the  storage  of 

-  the  arrays  COL  ft;  ]  ,  t  =  1,2,...  ,n  which  have  at  most 
m  components . 

-  the  n-arrays  N  ,  M  and  MI  ,  J 

-  the  "dynamic"  arrays  (at  most  m  components  each): 

1  =  1(0;  ] ,  1(1;  ] ,...,I[n-l;  ]  . 

-  the  current  linear  programming  tableaux  A  ,  which  is  at 
most  (n+m)  by  (rri-1)  . 

TREE:  1  -  Set  t  =  1  arid  J  =  0  ;  the  initial  tableau  stems 

from  (l);  1(0;  ]  =  1  ;  MI(0]  =  number  of  elements  in  1 

Mf  0]  =1 

-  Set  COL  [1]:  =  first  column  of  the  original  matrix  A 
(corresponding  to  the  first  non-basic  variable  x^) 

-  Set  N[l]:  =  j  =  ion-basic  index  ot  the  first  column 

-  Delete  the  first  column  from  the  matrix  A  ,  to 

obtain  eh*  ament  matrix  A  (which  is  thus  an  Ml[0J 
by  n  array) 

2  If  (n-t)  -  2  Lher  use  FACE2D  and  go  to  BACKTRACK. 

3  Apply  MINSET  to  the  current  tableau  A  ;  -»  i|t;  ] 

-  Set  MIft]  =  number  of  elements  in  the  set  lft;  ] 

-  Set  Ml  t]  -  0 

4  If  Mlft]  =  1  ,  apply  SIMFACE  and  go  to  BACKTRACK. 

5  If  M[t]  =MI[tl  ,  go  to  BACKTRACK  ; 
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6  -  Set  M[  t]  :  —  Mf  t]  +  1  ; 

7  Take  the  M[ t] th  component  j  of  the  array 
If  t;  ]  i.e. ,  j:  -  I[t;M[t]] ; 

8  Make  x^  non-basic  (if  it  is  not  already) 
preserving  primal  feasibility. 

When  pivoting  is  required,  one  must  transorm  also  the 
columns  stored  in  the  arrays  COLfs]  ,  for  all  s, 

1  <  s  <  t  ;  this  can  be  done  by  forming  a  full  tab¬ 
leau  TAB  =  [COLfl],  COLf 2] , . . . ,COL[ t] ,A}  ;  TAB 

is  then  transformed  by  pivoting  (pivot  in  A  )  and  the 
new  columns  1,...  are  stored  again  in  COLf 1] , . . . ,COL[ t] 

-  Set  t:  =  t  +  1  ; 

-  Set  N[  t]  :  =  j  ; 

9  -  Set  COLft]:  =  column  of  the  current  tableau  A  correspond¬ 

ing  to  the  non-basic  index  j 
-  Set  A  =  matrix  obtained  from  the  current  matrix  A  by 
deleting  the  column  j  ; 

10  Go  to  step  2. 

BACKTRACK: 

11  Set  t  =  t  -  1  . 

12  If  t  =  -1  then  STOP! 

13  Adjoin  the  column  COLft]  to  the  current  matrix  A  , 
thus  forming  a  new  (augmented)  matrix  A  . 


10. 


3.2  We  need  the  following  special  procedures  in  the  above  algorithm. 
Two-dimensional  faces:  FACE2D. 

Two-dimensional  faces  of  the  convex  polyhedron  P  can  be  treated 
separately,  because  a  straightforward  sequence  of  pivots  determines 
all  their  vertices  and  consequently  also  their  1 -dimensional  faces. 

FACE2D: 

1  Find  a  basic  feasible  solution,  i.e.,  two  non-basic  indices 

y 

-  Set  R  =  I [n-2 ; ]  and  J  =  0 

2  Register  the  corresponding  vertex  characterized  by  the  non- 
basic  index  set  VERTES  15  {N[1J,  N [2] , •  •  •  ,N [n— 2} ,iij,n2  }. 


3  Choose  a  non-basic  index  j  =  n^  or  n2  ; 

-  if  both  and  belong  to  J  then  STOP! 

4  Set  J:  -  J  v  {j}  and  R  =  R  -  {j} 

5  Find  the  basic  index  i  c  R  which  is  to  leave  the  basis 
when  the  non-basic  index  j  enters  the  basis  (in  order 
to  maintain  primal  feasibility). 
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Pivot  (on  a_  and  go  to 


2.) 


3.3  Slmplical  faces:  SIMFACE 

When  the  minimal  set  1  of  a  k-dimensional  face  F  consists  of 
(k+1)  elements,  F  is  a  simplex  and  its  (k+1)  vertices  can  be 
immediately  determined,  without  resorting  to  subfaces  of  F  . 

Lat  I  -  I  (i^  ,i2 , . . . ,  in— |t>  ( J]_  ’  ^2  ’  *  *  *^k+l^  * 


SIMFACE:  1  Find  a  first  basic  feasible  solution,  i.e. 

k  non-basic  indices  n^  , . . . ,n^  €  I  and 
store  VERTEX  =  {N[ 1] ,N[2] , . . . ,N[n-k] ,nx , . . . ,nRj 

2  For  each  i  =  1 . k  ,  generate  a  new  non-basic 

array  VERTEX  from  the  one  obtained  in  step  1 
above  by  replacing  the  index  by  the  last 

remaining  element  n^^  *  I  •  Clearly  n  new 
arrays  VERTEX  are  generated  in  this  fashion. 

The  procedure  SIMFACE  only  determines  the  vertices  of  F,  but  if  all 
substances  of  F  are  desired,  they  can  be  obtained  in  a  similar  manner. 

3.4  The  algorithm  TREE,  BACKTRACK,  FACE2D,  and  SIMFACE  determines  a 
non-redundant  list  of  non-basic  indes  sets  VERTEX  which  contains 
all  the  vertices  of  P. 

Proof :  The  procedure  TREE  generates  all  sets  of  indices 

N  *  £N[ 1] ,N[2) ,...,N[t]}  such  that  N[i]  e  I[i;  ]  , 
for  all  1  <  i  <  t  and  for  all  t,  0  <  t  <  n-2  ; 
for  all  t  >  0  ;  for  t  =*  n-2  the  procedure  FACE2D 
takes  over  and  finds  all  vertices  in  that  face  ; 
occasionally  for  0  <  t  <  n-2  ,  the  procedure  SIMFACE 
will  do  the  same.  Hence  all  basic  feasible  solutions 
of  P  with  non-basic  variables  in  the  minimal  set 
I  (of  P  )  are  generated  because  by  construction. 
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I  =  I[0;  ]  3  i[l;  ]  ^  . . .  o  l[t;  ]  always  holds  true. 

Furthermore  FACE2D  and  SIMFACE  determine  only  basic 
feasible  solutions. 

Moreover  this  list  is  non-redundant ,  because  step  6 
of  TREE  guarantees  that  the  same  index  set  N  is 
generated  once  and  only  once.  q.e.d. 

3.5  Degeneracy :  When  some  faces  (or  vertices)  of  P  are  degenerate » 
it  may  happen  that  geometrically  identical  faces  of  P 

are  algebraically  represented  by  different  sets  N  : 
such  faces  will  then  appear  more  than  once  in  the  arbor- 

escence  (but  each  time  with  a  different  index  set  N  , 
i.e.,  for  instance 

(x  I  x  c  F1]  =  {x  |  x  e  F 2) 

with  F^  =  F(i1>i2 , . . . jig) 

F2  —  ^ ( J j » ^2  *  * ' *  * Jq ^ 
an<^  {.  1^ » ’  *  *  *  ’  ^  { j  j » ^2  *  •  •  •  > 

This  redundany  is  due  to  degeneracy  in  the  arbor- 
escence,  i.e.,  to  the  fact  that  some  faces  F(i,,...,i  ) 

1  q 

are  "over-determined"  in  the  sense  that  3i  i  {i^,...,i^} 
such  that  F(i^,...,i^)  C  {x  /  x^  =  0}. 

Such  an  index  i  is  identified  by  the  procedure 
MIN  SET  and  one  has 

i  t  I(it , . . . ,i  )  but  i  e  I  (i....,i  ). 
i  q  o  1  q 


Hence  all  the  characterizations  by  different  index  sets 
of  the  same  face  F  can  easily  be  obtained 
from  the  index  set  D  *  [i  ,...,i  }  u  I  (i,)  o  I  (i.jOu... 

r  1  q  OL  O  i  4 

. . .  ul  (i_,...,i  );  One  has  by  definition 
o  l  q  J 

Dp  =  {  i  j  xt  =  0,  Vx  e  F} C  M 

and  Aff(F)  =*  [x  I  x.  =  0,  Vi  e  D  }.  The  corresponding 
linear  system 

x.  =  a.  +  £  a  x  =  0  ,  i  e  D.C  M 

1  \>  jcN  J  F 

has  rank  q  (<  n) ;  it  is  over-determined  in  the  sense  that 

contains  d  >  q  linearly  dependent  equations  (d  =  number 
r  F 

of  elements  in  D._)  whenever  one  of  the  sets  I  /  0. 

F  o 

Every  set  of  q  linearly  independent  rows  represents 
a  characterization  of  F  which  appears  in  the  arborescence 
and  causes  redundancy;  however,  because  the  main  subroutine 
TREE  generates  faces  in  a  lexicographically  increasing 
order,  redundancy  can  easily  be  eliminated  by  the  following 
additional  bookkeeping: 

Dl  -  Generate  the  first  (lexicographically)  degenerate 
representation,  i.e.,  the  next  redundant  node  to  be 
encountered  by  the  procedure  TREE; 

D2  -  Store  the  lexicographic  order  rank  r  of  that  node; 

D3  -  When  BACKTRACK  attains  the  order  r  (terminate  that 
redundant  branch  and)  go  to  Dl  in  order  to  update  r; 
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Note  that  the  degeneracy  sets  ere  generated  in  increasing 
lexicographic  order;  and  the  step  Dl  does  not  require  any 
additional  computations;  moreover  step  3  does  eliminate  redun¬ 
dancy  because  the  degeneracy  sets  newly  discovered  by  the 
procedure  TREE  all  belong  to  higher  orders  than  the  current 
one  (i.e.,  one  never  finds  out  in  a  later  phase  that  some 
previously  enumerated  faces  were  redundant). 
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Section  4:  Some  applications  of  a  Facial  Optimization  Method  (FOM)  for 
linearly  constrained  mathematical  programing. 

The  Facial  Optimization  Method  (FOM)  is  of  the  branch  and  bound  type, 
where  the  branching  procedure  Is  generated  according  to  the  facial  structure 
of  the  constraint  polyhedron  P  ;  in  this  manner  primal  feasibility  is 
always  present  in  all  the  nodes.  It  should  be  noted  that  the  arborescence 
of  section  2  still  contains  some  useful  flexibility  in  the  choice  of  the 
order  in  which  the  variables  are  to  be  "blocked"  at  value  0  (as  non¬ 
basics);  thus  the  usual  preference  rules  of  branch  and  bound  algorithms  also 
apply  here.  Since  in  every  way  the  facial  arborescence  resembles  ocher 
branching  trees,  the  usual  bound  estimators  al6o  apply;  clearly  the  effici¬ 
ency  of  the  facial  approach  as  compared  to  other  branch  and  bound  algorithms 
depends  on  P  and  its  structure. 

In  a  vertex  to  vertex  optimization  method,  where  additional  properties 
must  be  checked  as  in  4.1,  one  must  be  able  to  identify  all 
neighbouring  vertices  of  a  given  vertex  x  ,  i.e.  find  all  extreme  rays 
of  a  (degenerate)  cone.  Consider  a  basis  at  x 

x  =  x  +  I  a  x  >  0  ,  Vi  '  M  (4) 

1  1  jcN  J 

with  non-basic  set  M  .  For  the  rays  stemming  from  x  ,  only  the 

planes  that  are  tight  at  x  are  of  interest,  i.e. 

M’  =  (l  }  i  e  N  ,  or  xt  =  0)  C  M 

The  problem  of  finding  the  extrema  rays  of  the  cone  C 
C  -  {x  |  x.  >  0  ,  Vi  e  M’}  is  equivalent  to  that  of  finding  all  extreme 


V  X. 
jeN  J 

corresponding  system  (1)  can  be  obtained  in  one  pivot  from  (4). 
a  search  method  of  this  type  is  developed  for  concave  programs. 


points  of  the  (n-1) -dimensional  polyhedron  P  *  (x  e  C 


=  1};  the 
In  [3] 


4.1  Zero-one  progranming:  Here  the  polyhedron  P  is  a  subset  of  the  unit 
cube;  hence  every  feasible  integer  solution  is  a  vertex  of  P  . 

Starting  at  the  L.P.  optimal  solution  x  ,  a  branching  process  which 
determines  (not  necessarily  explicitly)  all  neighbours  of  x  ,  then  all 
neighbours  of  these  neighbours,  etc.... will  terminate  when  an  integer 
vertex  has  been  found  and  all  neighbours  generated  so  far  furnish  a 
lesser  value  of  the  objective  function.  Depending  on  the  structure  of 
P  ,  this  branch  and  bound  algorithm  can  use  efficiently  the  facial 
structure  of  F  ,  particularly  for  polyhedra  P  where  x  always  has 
an  integer  neighbour. 


4.2  General  quadratic  programming:  Suppose  we  want  to  optimize  (mximize 

or  minimize)  a  quadratic  objective  function  (not  necessarily  convex  or 
concave)  on  a  polyhedral  set  P  .  It  can  be  shown  that  the  candidates 
for  the  optimum  are  either  vertices  of  P  or  (more  generally)  optimal 
solutions  of  the  convex  quadratic  problems 

minimize  f(x)  ,  s.t.  x  e  F 


where  F  is  a  k-dimensional  face  of  P  such  that  the  function  f 
remains  convex  on  F  . 

The  facial  decomposition  method  leads  in  this  case  to  an  efficient 
algorithm  for  general  quadratic,  linearly  constrained  problems  which 
is  presented  in  (4  ].  ,, 

Of  course  FOM  can  be  applied,  conceivably,  to  many  other  special 
optimization  problems  over  polyhedral  sets  (such  as  concave  programming. 


general  non-ltnear  programming,  etc....)  but  it  is  not  clear  at  this 
point  where  and  when  it  furnishes  an  efficient  approach  in  such 
general  context. 

Conclusions 

The  algorithm  for  generating  all  vertices  of  a  polyhedron  P 
presented  in  this  note  develops  a  facial  arborescence  which  may  be  con 
sidered  minimal,  since  every  node  is  generated  exactly  once,  and  corre 
sponds  to  a  feasible  face  on  the  boundary  of  P  .  In  particular,  one 
can  generate  in  this  arborescence  the  set  of  vertices  of  P  (i.e. 
feasible  basic  solutions)  without  generating  any  other  basic  solution, 
as  is  the  case  for  methods  [1]  and  [2].  One  may  argue  that  the 
additional  pivoting,  necessary  for  the  construction  of  the  minimal 
sets  outweights  this  theoretical  property;  in  view  of  the  exponential 
growth  of  the  face  structure,  however,  this  should  not  be  expected 
to  hold  true  in  the  general  case.  A  program  requiring  a  modest 
amount  of  storage  (i.e.,  where  practically  only  space  for  A  is 
needed)  is  being  implemented.  When  degeneracy  occurs,  a  subtree 
of  the  facial  arborescence  can  be  eliminated  by  additional  bookkeeping 
The  remaining  arborescence  is  in  one  to  one  correspondence  with  the 
k-dimensional  faces  (n  >  k  >  0)  of  P  for  all  k.  Precise  comparisons 
in  this  context  are  not  easily  formulated  mathematically  and  experi¬ 
mentation  is  currently  under  way  with  randomly  generated,  non- 
structured  matrices  A  of  a  size  up  to  20  by  10. 
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